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Abstract 

The parametric excitation of pairs of electron surface waves (ESW) in the 
interaction of an ultrashort, intense laser pulse with an overdense plasma 
is discussed using an analytical model. The plasma has a simple step-like 
density profile. The ESWs can be excited either by the electric or by the 
magnetic part of the Lorentz force exerted by the laser and, correspondingly, 
have frequencies around oj/2 or to, where u is the laser frequency. 
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I. INTRODUCTION 



Surface wave (SW) excitation has been widely studied in the past as a mechanism for 
electromagnetic (EM) wave absorption in highly inhomogeneous plasmas in several regimes 
and for different target geometries |L],^| and, more recently, in the specific case of solid 
targets irradiated by intense, ultrashort laser pulses H,|J. The key problem in the linear 
mode conversion of the EM wave (the laser pulse) into an SW of the same frequency u is 
that the phase matching between the two waves is not possible at a simple plasma-vacuum 
interface. The reason for this is that for an SW vave the wavevector along the plasma surface 
k s is larger than u/c, and thus cannot be equal to the component of the EM wavevector in 
the same direction k t = (oo/c) sin 9, where 9 is the incidence angle of the laser pulse. Linear 
mode conversion of the laser pulse into an SW is possible in specifically tailored density 
profiles, e.g. for a "double-step" density profile (i.e. at the interface between an underdense 
plasma region where uj > u p and an overdense region where uj < uj p , with uj p the plasma 
frequency), or if a periodic surface modulation with wavevector k p = k t — k s exists, such 
that the matching condition for the wavevectors is 

k t = k s + k p . (1) 

In practical applications, corrugated targets are used to optimize SW generation at a given 
incidence angle. Other possible schemes of SW excitation are discussed, e.g., in Refs. 

In this paper, we discuss a different mechanism of electron surface wave (ESW) excitation 
by ultrashort laser pulses, based on the generation of two counterpropagating ESWs. The 
basic idea is as follows. An intense laser pulse impinging on an overdense, step-boundary 
plasma drives electron oscillations at the frequencies uj and 2uo by the electric and the 
magnetic part of the Lorentz force , respectively. The electron response may be viewed as 
a superposition of "forced" surface modes with frequency uj = uj or 2o>, respectively, and 
wavevector k$ = (u)q/c) sin 9. Each of these forced modes may "decay" parametrically into 
two ESWs (labeled "+" and "— " , respectively) provided that the matching conditions for a 
three-wave process hold: 
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k = k + + k_, (2) 
ujq = uj + + uj-. (3) 

No pre-imposed target modulation is required to satisfy these relations, as depicted in Fig.[TJ. 
The ESW frequencies may be written as uj± = ujq/2 ± 5uj. Thus, the "decay" of the mode 
excited by the electric (magnetic) force leads to the generations of ESWs with frequencies 
around uj/2 (uj). We call this process "two-surface wave decay" (TSWD) in analogy to the 
well-known "two-plasmon decay" parametric instability in laser-plasma interactions. 

Earlier investigations of parametric and three-wave nonlinear processes involving ESWs 
have been previously carried out in different regimes ||||. In particular, the parametric 
excitation of pairs of subharmonic ESWs in a magnetized plasma by an external EM wave 
was studied in Refs. [[/J. These studies were restricted to the electrostatic limit in which 
the ESW frequency is uj s ~ uj p /\/2 (see || and eq.([TT|) below) and does not depend on 
the wavevector k s in a semi-infinite, step boundary plasma (e.g. with a density profile 
rii = riiO(x), where Q(x) is the Heaviside step function). In the electrostatic limit the phase 
matching with an impinging EM wave is possible only for particular target geometries where 
the frequency of electrostatic SWs depends on k s (e.g. a plasma slab of thickness d where uj s 
depends on the product k s d fBJ). In our model, we assume a semi-infinite geometry but we do 
not use the electrostatic approximation, and therefore in general the frequency uj s of ESWs is 
significantly lower than uj v I \/2 and depends on k s , approaching the EM limit uj s ~ k s c when 
uj s — > 0. This allows the matching conditions (§) to be satisfied in a semi-infinite plasma. 
It is interesting to notice that in Ref. it has been shown that two counterpropagating 
ESWs generate nonlinear radiation at a frequency that is the sum of the ESW frequencies; 
this process may be regarded, in some sense, as the inverse process of the TSWD. 

One of the main results of this paper is that the intense v x B force at 2uj of the 
laser pulse may drive a 11 2uj — > uj + uj" decay with a growth time of a few laser cycles. 
In this paper, we focus on this case because of its direct relevance to the interpretation 
of recent numerical simulations in two dimensions (2D) for normal laser incidence, which 
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have provided evidence for this process [10|. Note that the surface oscillations observed in 



Ref. [ECO] and theoretically studied in the present paper are different from the grating-like 



surface inhomogeneities induced by the magnetic force that were studied in Ref. JO]; these 
latter "non-parametric" surface structures oscillate at the same frequency 2u of the driving 
force and are generated at oblique incidence only. In the case of a p-polarized laser pulse at 
oblique incidence, we show below that the parametric process may be driven by the electric 
force leading to the generation of subharmonic ESWs around the frequency u/2. 

The primary aim of this paper is to describe the "basic principle" of the TSWD as 
a possible route to the generation of surface waves and to outline an analytical model 
for the TSWD. Therefore, for the sake of simplicity, we use a cold fluid, non-relativistic 
plasma model that can be tackled analytically. However, we believe that the TSWD is a 
general process that may occur for laser and plasma parameters beyond the limits of the 
approximations used in the present paper. Relativistic and kinetic effects, the inclusion 
of damping mechanisms and the fully nonlinear evolution of the TSWD are left for future 
investigations. 

II. THE MODEL 

In our model, we consider a cold plasma with immobile ions and a step-like density 
profile Tii(x) = riiQ(x). The laser pulse has frequency u, linear polarization and wavevector 
k = (uj/c) (x cos 9 + y sin#). Translational invariance along z is assumed. For both s- and 
p-polarizations, all fields oscillating at the frequency u can be determined inside and outside 
the plasma from general Fresnel formulae for refraction and transmission at a boundary 
12|| between vacuum and a medium with an (imaginary) refractive index n given by n 2 = 



1 — u 2 /u 2 < 0, Up = (Annie 2 / me) 1 1 2 being the plasma frequency. 

In what follows we study the excitation of "H" ESWs with the magnetic field in the 
z direction. Thus, we deal with an effectively 2D geometry, in which B = Bz and the 
Maxwell-Euler equations read 
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V-E = 47re(n -n e ), (4) 

VxE = --d t B z , (5) 
c 

47T 1 

V x B = — j + -dtE, (6) 

c c 

m e d t v = — m e v ■ Vv — e + — x (7) 

and the current density is given by j = — en e v. 
We adopt the following expansion for all fields 

F(x,y,t) = F { (x) + eF {x,t-y sin 6/c) + e 2 [/+ (x, y, t) + /_ (x, y , £)] , (8) 

where e is a small expansion parameter and / stands either for the electron density or the 
velocity or for the EM fields in the (x, y) plane. In this expansion, Fi represents unperturbed 
fields of zero order (e.g., ni). The term F represents the "pump" field at the frequency u , 
that is taken to be of order e and is written as 

In what follows the cases with uj =uj and Uq = 2u are studied separately. The last term in 
(§) is the sum of two counterpropagating surface modes, which are assumed to be of order 
e 2 and are written as 

f±= l -[ft ± \x)e^y- i ^ t ]+c. C .. (10) 

Using expansion (||), the coupling between the pump and the surface modes is of order e 3 . 
Thus, from eqs.(f|-(5|-|6|-|7]) we obtain to order e the pump fields. To order e 2 , we obtain a 
set of linearized, homogeneous 2D Maxwell-Euler equations which have ESWs as normal 
modes. Note that terms such as v • Vv in the Euler equation are also of order e 2 , but are 
non resonant with the ESWs at the frequencies uj±. Such terms represent a source term for 
harmonic oscillations at 2uj and can be neglected in our model. 

ESWs are "H" waves which can propagate in a dielectric medium along a layer of dis- 
continuity of the dielectric function, with the latter changing sign across the boundary [13 . 



For a cold, step-boundary plasma there is no volume charge density perturbation associated 
with SWs, i.e. V • E = — AneSrig = 0. At a vacuum-plasma interface, electrons do not enter 
the vacuum side, but are "stopped" at the interface (x = 0) forming a surface charge layer. 
The dispersion relation for an SW of frequency u s and wavevector k s is 



u2 = uj 2 - uj = u 2 a s - 1 



(11) 



We have set for convenience a s = uj 2 /uj 2 . The dispersion relation ( |TT1) is shown in Fig.|T]. 
The EM field envelopes (see eq. fllOD ) are given by 



E s y (x) = E s y (0) [Q(-x)e q < x + e{x)e~ q > 3 
B s z {x) = l -^E s y (0) \®{-x)e q < x + e{x)e- q > 3 

E s x = ikE s y (0) Q(- x )—-e(x) 



-q > x 



q< 



q> 



(12) 



where g> = (u> s /c)(a s — lj\/a s — 2), q < = (u s /c)(l/ \Ja s — 2). The velocity fields are found 
from v s = (ie / 'm e uj S )E S . 

All feedback effects of the nonlinear coupling on the pump fields are neglected. The 
coupling terms that lead to the excitation of the ESW with frequency u± may be represented 
by the nonlinear force 



t+ L) (x)e ik±y - w±t 



m e (v T • Vv + v • Vv T ) + -(v x B T + v x B ; 



(13) 



In the term in square brackets on the r.h.s of ( |T3| ) only "resonant" terms with the same phase 

of 4 iVL) 

are included. Obviously, these terms exist if the matching conditions ([$]) hold. 
The nonlinear coupling leads to the parametric excitation and growth of ESWs. We thus 
let the ESW envelopes vary slowly in time, i.e. 



f±{x) -> f±(x,et) 



(14) 



To evaluate the growth rate, we use an energy approach. The surface energy density of the 
ESW that will be needed in the calculations is given by 
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Z7T J-iv/k J-oo x ' 

where the brackets denote average over one laser cycle and u± and are the volume 
densities of the kinetic and EM fields energies, respectively: 

^ = ^N', (16) 
M « = ±(|E ± | 2 + |B ± J 2 ). (17) 

The integral over ?/ yields a factor 1/2. The kinetic energy contribution vanishes for x < 0. 



Using eqs.(12) and v±, x (0 + ) = f±, y (0) /^a — 1, and finally summing the two contributions 
we obtain for the total surface energy U± = M±\v± >y \ 2 /2 where 

m e nica±(a± - 2) 1 / 2 (a| - 2a± + 2) 
M± = (^1? ' (18) 

To obtain an expression for the temporal variation of U± due to the nonlinear coupling 
forces, we write the Euler and the Ampere equation for the ESW fields keeping terms up to 
order e 3 and disregarding terms that are non-resonant with the oscillation: 

m e d t v ± = -eV ± + e^ L \ (19) 

V x B ± = — (-en oV± + eJ ( ± L) ) + -d t E ± , (20) 
c K ' c 

where is given by eq.(|l3|), and J± Vi ' ) = — e5n < f°^v±. Using eqs. (|l9| - p0|) together with 
Maxwell's equations we obtain for the temporal variation of the energy densities 

d t uf ] = -en D v ± • E± + eff L) • v±, (21) 
dtu^ = -V • S± - ejf L) ■ E±, (22) 

where S± = (c/4ir)'E± x B± is the Poynting vector of the ESW. We integrate the sum of 
the energy densities over space and average over time. Noting that the total flux of S± 
vanishes because of the evanescence of the ESW fields, we finally obtain the equation for 
the evolution of the total energy U of the ESW: 
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d t U± = — dy dxdt{——\v±\ + — ( E± + B± 

J-Tr/fc JO \ 2 87T V 



J-n/k 

9ir /"+oo r+n/k , , . . 

/ rfx / dy ( v± . ( e ^)E ± + ff-)) . (23) 



The integral extends for x > only since there are no electrons for x < 0. 

From eq . (p3|) , performing the integral and neglecting high order terms one obtains two 
coupled equations in the general form 

^7^dt\v±\ 2 = A u\v + \\v + \S± sin if, (24) 

where A is the amplitude of the pump mode in dimensionless units, S± are positive factors 
depending on u p /u and 9, and is obtained from the phase factors of the pump and ESW 
modes as <p = (p + + <p- — <fio- From (E3) one easily obtains 



8^\v ± \ = ±^ V \v ± \. (25) 

Setting sin (p — 1 simply corresponds to finding the relative phase of the growing modes with 
respect to the phase of the pump mode. The growth rate of ESWs is then 



S+S. 



T = A ujJ—± . (26) 

In the next section we use the analytical method outlined above to study TSWD in two 
cases of particular relevance: TSWD driven by the v x B force for normal laser incidence 
(which we name "2cj — > uj + u n decay), and TSWD driven by the electric force for oblique 
laser incidence and p-polarization ("a; — > u/2 + ui/2" decay). 

III. "2u -» u + w" DECAY 

To discuss the u 2u — > u> + u" decay, we restrict ourselves for simplicity to the case of 
normal incidence, as was done in the simulations reported in Ref. |[L0|| . Thus, the ESW 
frequency uj± = uj and k + = —k- The case of oblique incidence is indeed very similar, the 
main difference being that the matching condition k + + k_ = (2u/c) sin# holds, causing a 
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shift Sou 7^ between the two ESWs. It is found that the growth rate decreases monotonically 
with increasing incidence angle, since the pump force has a maximum at normal incidence 

0- 

The laser wave can be represented by a single component of the vector potential, A z = 
A z (x, t), that for x > is given by 

A z (x,t) = A z (x)cosut = A z (0)e- x/la coswt, (27) 

where l s = c/(u 2 - u 2 ) 1 ' 2 is the screening length. Imposing boundary conditions for the 
incident and reflected waves one finds 4,(0) = 2Ai(ujl s /c)(l + u 2 l 2 Jc 2 )~ 1 / 2 , where A { is the 
amplitude of the incident field. 

Electrons perform their quiver motion in the z direction. Thus, the magnetic force term 
is in the x direction and has a secular term (Ooj), named the ponderomotive force, and an 
oscillating term (2u) that, in what follows, we simply name the v x B force. The secular 
term corresponds to radiation pressure and creates a surface polarization of the plasma. 
The v x B force drives a longitudinal, electrostatic oscillation that acts as a pump mode 
for TSWD. In the expansion (^), the pump amplitude is supposed to be of order e. As will 
be shown below, coupling between ID and 2D fields occurs only in the overdense plasma 
(x > 0). Thus, since the v x B force is quadratic in the laser field, the expansion (H) 
also implies a(0) ~ e 1 / 2 , where a(0) is the (dimensionless) laser amplitude at the surface 
of the plasma, e.g a(0) = (eA z (0)/mc 2 ). For overdense plasmas a(0) ~ (u>/uj p )ai < a i: 
where a, = (eAi/mc 2 ). This yields e ~ [uj/u p ) 2 a 2 = (n c /n e )a 2 , where n c = m e uj 2 /Aiie 
is the "critical" density. We therefore expect our expansion procedure to be valid even at 
relativistic fields amplitudes a« > 1 for high enough plasma densities. 

We now derive the solutions for the ID pump fields at 2uj. The electron oscillation 
velocity is obtained from the conservation of canonical momentum as m e v z = eA z /c. For 
normal laser incidence, the longitudinal v x B force is given by 

2 

- -v z B y = e —A z d x A z = F°{x){l + cos2cjt), (28) 

c m e c z 
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where we have set F°(x) = {rrieC 2 /2l s )a 2 e 2x l ls = F°e 2x l la . To order e, one obtains the 
following equations for the longitudinal, electrostatic motion : 

rriedtV^ = -e{Ef ] + E {2w) ) + F°(x)(l + cos2cut), (29) 
dM 0) + 4 2UJ) ) = -4*e(6n<p + 8nf^), (30) 
d t 5nM = -nod x V^\ (31) 

All fields in the equations above decay inside the plasma as exp(— 2x/l s ). The secular part 
simply gives eE^(x) = F°(x) and 5n^ = —d x F°(x)/4:7re 2 = F°(x)/2ire 2 l s . For the motion 
at 2a; one obtains 

T>(2^) _ ~ tF ° ( V_ \ p -2x/l s ( on\ 

K -com e D[co 2 ) e ' {62) 
6h^= ni J° e~ 2x / l % (33) 

eE^ = ^-e~ 2x/h . (34) 
2D 

The denominator D is given by 

D=l -% = 1 -^, (35) 

which shows the well-known resonance at n e = 4n c due to excitation of plasmons with 
frequency u p = 2uj by the v x B force. 

The difference between the total number of electrons and ions for x > is, to order e, 

F° / cos2wr 



AA/" (l>0) 



/> K »» + , n f)^( 1 + ^). (36) 



The fact that AN^ X>0 ' > during most of the oscillation implies that, due to compression 
from the ponderomotive and v x B forces, electrons leave behind a charge depletion layer 
of thickness £ = AN^ X>0 ^ / 'rij. We note that £ is of order e, thus to lowest order it is correct 
to treat the charge depletion layer as a surface layer. Heuristically, the oscillating behavior 
of £ describes the plasma "moving mirror" effect ||15|| : due to charge depletion the laser 
is reflected at x = ( rather than exactly at x = 0. This leads to the appearance of high 
harmonics in the reflected light. 
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Since \D\ < 1, there is always a time interval in which electrons are pulled into vacuum 
forming a cloud of negative charge. This interval is very short for n e /n c ^> 4, i.e. D ~ 1, 
but for lower densities it is relevant in the interaction process. The motion in vacuum 
is strongly anharmonic and more difficult to solve analitically than the motion inside the 
plasma. However, it is found in section |TTT] that the surface modes gain energy during the 
phase of electron motion inside the plasma only, so that the expressions of fields for x < 
are not needed. Nevertheless, we have to assume that the density of the electron cloud for 
x < is low enough for the laser propagation not to be affected. 

Using eq.(32), the NL force ( |i~3"D may be written in terms of the oscillation velocity V^ 2 ^ 
as follows: 



cNL 



m e (y^\x)d x v T + v x ^d x V^\x)*) - -V^(x)B ZjT y 



(37) 



Because of the symmetry with respect to inversion of the y axis, the two ESWs may 
differ only by a phase factor and have the same amplitude and energy. Inserting the force 
( P7|) in the energy equation (^) and performing the integral in y we obtain for the ESW 
energy U = U± 



d t U = \ E r°° dx Ki ■ (e$n e ^ (x)E* , - m e n v* x .^V^ (x) 



l=+k,-k' 

-m e n V x ^ ) (x)d x r i + n -V x ^\x)B* z _,y ) + c.c. 

c ' 



(38) 



We rewrite all fields as functions of V^ 2 ^ and v±. We obtain, after some algebra 



m e rn 



-m. I vS 2u3) 



2{l + q> l 

All the terms of eq.(|39D are proportional to 



1 + q> l s )v* + v*_ + 2v*v* 



(39) 



3? (y( 2uj) v* + v*_) = \V x {2uj) v 2 ± \ cos(0 + tt/2), 



(40) 



where = (f> + + 0—, <j)± being the phase angles of v±, and we used eq.(^2]) and the fact that 
F° is real and positive. Thus, the phase of the growing ESWs with respect to the v x B force 
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is such that <fi = — 7r/2, the value for which the growth rate is positive and has a maximum. 
The overlap of the two ESWs gives 

^+ iky—iu)t I ^~ —iky—iujt l^i I —iujt f iky+i<f>+ _i_ —iky+i<f>-\ 

2 2 ~ 2 ^ / 

= |C ± |e-^- i7r/4 cos (ky + A0) , (41) 

where A<p = (<ft + — 0_)/2. This is a standing wave which has a temporal phase shift — 7r/4 
with respect to the v x B force of eq. (|28|) . The phase is such that, at a given position in y, 
the temporal maxima of V^ 2 ^ and v ^ overlap once for laser cycle. The angle A0 gives the 
location of the maxima of the standing wave on the y axis, which depends on the arbitrary 
choice of the initial phase. 

Using eq. ( ^0|) , we can rewrite the energy variation as 



d t U = m eni \V^ 



in) I 



l 2 + 



i + q>l 3 q>l s (i + q>h 



(42) 



Note that d t U > since g>/ s = y (a — l)(a — 2) _1 > 1, where a = oj^/uj 2 = n e /n c . 
Eliminating V^ 2 ^ as a function of a, and a, after some algebra the growth rate T is obtained: 



2 



[a 



- If' 2 



o;|a-4|[(a-l) 2 + l](a-2) 1 / 



2 



2(a- 1) 



a(l + q>l s ) ag>/ s (l + g>/ s 



- 4wa, 2 — 777-^ -ryr. (43) 

4 a|a-4|[(a-l) 2 + l](a-2)V2 v 7 

The leading contribution was highlighted in the last equality. 

A plot of the growth rate as a function of a is given in Fig.[| For moderately overdense 
plasmas, T is a considerable fraction of the laser frequency. The resonance at n e = 4n c 
is due to the resonant behavior of the pump field. The growth rate diverges also in the 
electrostatic limit when n e —>■ 2n c . However, in this limit the value of k tends to infinity, 
i.e. the wavelength becomes very small and one expects that this second resonance might 
be damped by thermal effects, which are neglected in our model and are left for future 
inevstigations. 

The comparison of spatial and temporal scales predicted from this model with the sim- 
ulation results shows reasonable agreement for laser and plasma parameters such that our 
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expansion procedure is valid |T0| . We note that simulation results for high intensities in the 
relativistic regime suggest that a 2uj —>■ uj + uj TSWD-like process still occurs and produces 
strong rippling of the plasma surface; the spatial scales are different from those predicted 
by our analytical, cold fluid model, but still of the same order of magnitude. 

We conclude this section by noting that, although we took the laser polarization along 
z, the whole derivation is valid at normal incidence for any polarization direction. The only 
difference in 2D geometry is that, for the laser polarization along y, the quiver oscillations 
at uj and the ESW oscillations at the same frequency overlap. This is observed in numerical 



simulations [14 



IV. "u ^u/2 + u/2" DECAY 

A case of particular interest is that of oblique incidence and p-polarization of the laser 
pulse, since these conditions lead in general to a stronger plasma coupling. In the preceding 
section we have found that the 2uj — > uj+uj decay is driven by the longitudinal component of 
the velocity field, as shown by eq. (p7|) . For p-polarization and 8^0 the electric field drives a 
strong longitudinal oscillation at the frequency uj; thus, one expects that the TSWD process 
leads to the generation of subharmonic ESWs, having frequencies around uj/2, with peak 
efficiency at some angle of incidence. As shown below, this picture is substantially correct. 
However, there is a small but non- vanishing growth rate for the uj — > uj/2+uj/2 process even 
at normal incidence. In fact, the pump force (|i~3| ) is now given by 

ff L) = -m e (v^d x VM + V^d x v T ) - me^VM + V y ^d y v T ) 

-^(V (w) x Bzp + Vzp x B (w) ), (44) 

and does not vanish for 8 = 0. 

The case of oblique incidence and p-polarization of the laser pulse may be tackled with 
an approach analogous to that of section |TT1[ The main physical difference with the previous 
case is that now the "pump" is a divergence-free velocity field at the laser frequency u, and 
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its amplitude is now proportional to the laser field rather than to the laser intensity. In 
applying the expansion procedure (§) again, now e ~ {uj/uj p )ai (appendix [A]). Thus, the 
limits of validity of this expansion are more restrictive in the present case: the expansion 
tends to be valid only for rather high densities or rather low fields. One must also note 
that temperature effects, neglected in the cold plasma approximation, are important for low 
fields. 

The calculation of the growth rate for subharmonic ESWs can be performed by the same 
method as in the preceding sections in a straightforward way. The details are reported in 
appendix |A[ As a final result, the growth rate V = T(a, 9) may be written as 



where Tb is the Fresnel factor for the magnetic field and the factor K(a, 9) scales weakly 
with density. A plot of T as a function of 9 for different values of a is reported in Fig||. 
The frequency shift 5uj of the two ESWs can be calculated as a function of 9 and a from 



Fig.|3], one finds that 5u ~ O.lu in conditions favorable to the uo — > uo/2 + uj/2 TSWD. 



The TSWD concept has been investigated analytically in two cases that clarify its basic 
features in the context of the cold fluid plasma approximation. This process may lead to the 
excitation of ESWs by an ultrashort laser pulse in a solid target without any special (e.g., 
grating-like) structure, and even for normal incidence and s-polarization. 

In both the cases discussed in this paper, we found a strong decrease of the TSWD growth 
rate for increasing densities. This does not necessarily imply that TSWD is not relevant to 
short pulse interaction with solid targets which have n e /n c ^> 1. In fact, important processes 
such as high harmonic generation or fast electron production are more efficient when the laser 
pulse interacts with a moderate density plasma. This is the case for most of the experiments 




(45) 



the matching condition k + + k_ = (uj/c) sin0. The result is shown in Fig.[|. Comparing with 



V. DISCUSSION 
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of short pulse interaction with solid targets, since a moderate density "shelf" is usually 
produced at the time of peak pulse intensity by target ablation and plasma hydrodynamic 
expansion during the leading edge of the short pulse or during the long prepulse. For 
instance, the ll 2u — ► uj + uj" decay appears to be very efficient exactly in conditions favorable 
for high harmonics production via the "moving mirror" effect [|T3|,|TBJ. Simulations fTU| 
suggest that TSWD in strongly nonlinear and relativistic regimes (beyond the limits of the 
analytical approach of this paper) may produce surface perturbations acting as a "seed" for, 
e. g. , electron instabilities leading to current filamentation |l7|Jl8| or detrimental distortions 
of the moving mirrors, and for Rayleigh- Taylor-like (RT) hydrodynamic instabilities occuring 
on the time scale of ion motion. In Ref. |H| it is estimated that the growth rate of the RT 
instability is much slower than that of TSWD. A different hydrodynamic instability with a 
surprisingly high growth rate has been reported in Ref. fL9f). 

The 2uo — > uo + uo and the uo — > uo /2 + uo /2 processes have been considered independently. 
This is appropriate since these are resonant processes which do not interfere with each other. 
In principle, both decays may occur during the interaction of an intense laser pulse with 
an overdense plasma. According to our analysis the 2uo — > uo + uo process appears to have 
a stronger growth rate. However, in the case of the uo —> uo/2 + u/2 process our analytical 
approach is valid for a very narrow range of parameters only. Either an extension of the 
present analytical approach or numerical simulations are needed to investigate the TSWD 
for a very intense, p-polarized laser pulse at oblique incidence. 



VI. CONCLUSIONS 

We have discussed the Two-Surface Wave Decay as a novel mechanism for the excitation 
of electron surface waves in the interaction of ultrashort, intense laser pulses with overdense 
plasmas. TSWD is based on the parametric excitation of a pair of ESWs. The "pump" force 
may either be the magnetic or the electric force of the laser pulse, and leads respectively 
to the generation of two ESWs with frequencies close to the laser frequency of half of it. 
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An analytical model for TSWD in the cold fluid plasma, non-relativistic approximation has 
been developed. The model supports the interpretation of recent simulation results [TIJ and 



suggests that TSWD may be of relevance in certain regimes of laser interaction with solid 
targets. 
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APPENDIX A: GROWTH RATE OF THE uj -> uj/2 + u/2 DECAY 

We now give the detailed derivation of the growth rate for the u — > u/2 + u/2 TSWD of 
section [TV} The pump fields are easily found with the help of Fresnel formulae and Maxwell 
equations. For instance, the magnetic field is given by 

B^{x,t) = B ^\0 + )e ikys ' me - x/lp - lujt + c.c. , (Al) 

where the screening length is given by 

*P = - ( 1 - -2 cos2 e ) V2 = ~ , 1 2a > (A2) 
w p \ / uj Va - cos 2 9 

and the magnetic field at the surface is given by the Fresnel formula 

B z (0+) _ 2n 2 cos# 
B z>i Vn 2 — sin 2 9 + n 2 cos 9 

= 2(a ~ 1)C " Sg =F B {0), (A3) 

(a — 1) cos 9 — i\f a — cos 2 9 

where B z i is the incident field amplitude in vacuum. The electric field components are found 
by using V • E (uj) = and V x E H = i(u/c)B^ + c.c. 

The nonlinear force is given by eq. (p7|) . Evaluating the spatial derivatives as d x \^ = 
-V^/I p , dyV^ = iktV, d y v± = ik±v±, and d x v± = —q±v± where q± = q > (u±), and 
keeping resonant terms only we find 
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r(NL) _ 
J±,x — . 



1 



f: 



(NL) _ rne 
±,v ~ a 



12 J ff~\j gC 

1 ~, 
v 

i(k t - k T )V^v^ y 



( V {uj) Bt + B { 



z u t,v 



e -(q T + l/l v )x + QQ 



(A4) 



{q T +l/l v )x 



] B* + B^v* 

x T,z ' z =F,a; 



+ c.c. (A5) 

We rewrite all "pump" terms as a function of V^ 03 ' by using t/. = ik t l p V y ^\ B^ 3 ' = 
(m e c/el p )(k^lp — 1)V1 ■ in a second step we rewrite all SW fields as a function of v± iV by 
using v± )X = (ik±/q±)v±^ y and b± )Z = —m e uj±(a± — l)v± ty /(ecq±). We thus obtain 

,2 



k lo 1 I 

-± (A* + fc^Zp + fc T ) + ^(« T - 1) + y 



— a v y u T,y^±> x(i 



7{NL) _ me™. t 



=F,1T 



, , / 1 hk T 
hip \ q T ~ y 



-{q T + l/l p )x 



(a T -l))+k T 



a v y D ^,y l ^±,y (i 



The relations u>±(a±—l) / q± = u±c^/a± — 2, k±/q± = 1/ \Ja± — 1, k t l p = sin6 l (a— cos 2 
may be used to simplify the expressions for Q± >x and Q± >y . We thus obtain 



Q 



±,X 



V 7 ^ 



-\k t {\ + q T l p ) + k T ] + ~-^Ja T - 2 + 



1 



Q±, y = hip [q T - j- 



^ x la T -2) +k T . 



(A6) 
(A7) 
(A8) 
(A9) 

0)-l/2 

(A10) 
(All) 



\ l P V"t 

In the energy equation (p3|) , since (faA^ = the integrand, eliminating w± iX , is given by 



v ± ■ £') = ^V^vlA, (^Q ± „ + Q ± ,„) + ex. 



Q± 



~-Q±,x + Q±, 



We thus find 



d t U± = n. 



-Q±,x + Q±, 



simp. 



(A12) 
(A13) 

(A14) 



1 + (q+ + q-)l p \q± 

Here, <p = (<j> + + </>_ — $), where $ and 0± are the phase factors of V y ^ and v±, y , respectively 
It is worth rewriting these equations in a more compact form. Introducing eij = eBi/m e ujc 
and using Fresnel's formulas to eliminate V y we obtain 
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Tfl 

d t U± = ni^a i c\J r B {9)\\v+,y\\v-,y\suiipG±(a,0), (A15) 



8 

where we have posed G±(a, 9) = G (a, 9)g±(a, 9) and 



Va - cos 2 9 

°° ~ (a-l)[l + ( g+ + g _)Z p ]' (A16) 



9± — lp \ —Q±,x + Q±.y 



a+-l)(a_-l) c V"± " 1 V"± " 1 



+ktl p \q^l p - 1 - ^=%y " V a T " 2J + Mp- (A17) 

Using the equations above we obtain the following coupled equations in the form (|4]) for 
the field amplitudes 

,y\ — ai^\v+y\\v^y\\J-B\G± sin ip, (A18) 

where fi± = 4M±/(m e nic). Thus, the growing modes have relative phases such that simp = 
1. The growth rate is given by 



T = aMF B \G > = aiw\r B \ ( —) K(a, 9) (A19) 

\ \n e J 

which is the formula (f£5|). A plot of T is depicted in Fig|| 
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FIGURES 

FIG. 1. Dispersion relation of electron surface waves (thick lines) eq.(^TJ), and the matching 
conditions for TSWD, eq.(||). 

FIG. 2. The TSWD growth rate T (normalized to af ui) for the 2ui — > ui + ui process at normal 
incidence, eq. (|43|) , as a function of a = n e /n c . 

FIG. 3. The TSWD growth rate T (normalized to aiui) for the u —* uj/2 + uj/2 process at 
oblique incidence, eq. (fj5[) , as a function of 6 and a = n e /n c . 

FIG. 4. The frequency shift 5uj for the uj — > to/2 + uj/2 process at oblique incidence, as a 
function of 9 and a = n e /n c . The dashed line gives the shift 8oj max = (1/2) sin 9 that is obtained 
for a — > oo. 
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Abstract 

The parametric excitation of pairs of electron surface waves (ESW) in the 
interaction of an ultrashort, intense laser pulse with an overdense plasma 
is discussed using an analytical model. The plasma has a simple step-like 
density profile. The ESWs can be excited either by the electric or by the 
magnetic part of the Lorentz force exerted by the laser and, correspondingly, 
have frequencies around oj/2 or to, where u is the laser frequency. 
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I. INTRODUCTION 



Surface wave (SW) excitation has been widely studied in the past as a mechanism for 
electromagnetic (EM) wave absorption in highly inhomogeneous plasmas in several regimes 
and for different target geometries |L],^| and, more recently, in the specific case of solid 
targets irradiated by intense, ultrashort laser pulses H,|J. The key problem in the linear 
mode conversion of the EM wave (the laser pulse) into an SW of the same frequency u is 
that the phase matching between the two waves is not possible at a simple plasma-vacuum 
interface. The reason for this is that for an SW vave the wavevector along the plasma surface 
k s is larger than u/c, and thus cannot be equal to the component of the EM wavevector in 
the same direction k t = (oo/c) sin 9, where 9 is the incidence angle of the laser pulse. Linear 
mode conversion of the laser pulse into an SW is possible in specifically tailored density 
profiles, e.g. for a "double-step" density profile (i.e. at the interface between an underdense 
plasma region where uj > u p and an overdense region where uj < uj p , with uj p the plasma 
frequency), or if a periodic surface modulation with wavevector k p = k t — k s exists, such 
that the matching condition for the wavevectors is 

k t = k s + k p . (1) 

In practical applications, corrugated targets are used to optimize SW generation at a given 
incidence angle. Other possible schemes of SW excitation are discussed, e.g., in Refs. 

In this paper, we discuss a different mechanism of electron surface wave (ESW) excitation 
by ultrashort laser pulses, based on the generation of two counterpropagating ESWs. The 
basic idea is as follows. An intense laser pulse impinging on an overdense, step-boundary 
plasma drives electron oscillations at the frequencies uj and 2uo by the electric and the 
magnetic part of the Lorentz force , respectively. The electron response may be viewed as 
a superposition of "forced" surface modes with frequency uj = uj or 2o>, respectively, and 
wavevector k$ = (u)q/c) sin 9. Each of these forced modes may "decay" parametrically into 
two ESWs (labeled "+" and "— " , respectively) provided that the matching conditions for a 
three-wave process hold: 
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k = k + + k_, (2) 
ujq = uj + + uj-. (3) 

No pre-imposed target modulation is required to satisfy these relations, as depicted in Fig.[TJ. 
The ESW frequencies may be written as uj± = ujq/2 ± 5uj. Thus, the "decay" of the mode 
excited by the electric (magnetic) force leads to the generations of ESWs with frequencies 
around uj/2 (uj). We call this process "two-surface wave decay" (TSWD) in analogy to the 
well-known "two-plasmon decay" parametric instability in laser-plasma interactions. 

Earlier investigations of parametric and three-wave nonlinear processes involving ESWs 
have been previously carried out in different regimes ||||. In particular, the parametric 
excitation of pairs of subharmonic ESWs in a magnetized plasma by an external EM wave 
was studied in Refs. [[/J. These studies were restricted to the electrostatic limit in which 
the ESW frequency is uj s ~ uj p /\/2 (see || and eq.([TT|) below) and does not depend on 
the wavevector k s in a semi-infinite, step boundary plasma (e.g. with a density profile 
rii = riiO(x), where Q(x) is the Heaviside step function). In the electrostatic limit the phase 
matching with an impinging EM wave is possible only for particular target geometries where 
the frequency of electrostatic SWs depends on k s (e.g. a plasma slab of thickness d where uj s 
depends on the product k s d fBJ). In our model, we assume a semi-infinite geometry but we do 
not use the electrostatic approximation, and therefore in general the frequency uj s of ESWs is 
significantly lower than uj v I \/2 and depends on k s , approaching the EM limit uj s ~ k s c when 
uj s — > 0. This allows the matching conditions (§) to be satisfied in a semi-infinite plasma. 
It is interesting to notice that in Ref. it has been shown that two counterpropagating 
ESWs generate nonlinear radiation at a frequency that is the sum of the ESW frequencies; 
this process may be regarded, in some sense, as the inverse process of the TSWD. 

One of the main results of this paper is that the intense v x B force at 2uj of the 
laser pulse may drive a 11 2uj — > uj + uj" decay with a growth time of a few laser cycles. 
In this paper, we focus on this case because of its direct relevance to the interpretation 
of recent numerical simulations in two dimensions (2D) for normal laser incidence, which 
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have provided evidence for this process [10|. Note that the surface oscillations observed in 



Ref. [ECO] and theoretically studied in the present paper are different from the grating-like 



surface inhomogeneities induced by the magnetic force that were studied in Ref. JO]; these 
latter "non-parametric" surface structures oscillate at the same frequency 2u of the driving 
force and are generated at oblique incidence only. In the case of a p-polarized laser pulse at 
oblique incidence, we show below that the parametric process may be driven by the electric 
force leading to the generation of subharmonic ESWs around the frequency u/2. 

The primary aim of this paper is to describe the "basic principle" of the TSWD as 
a possible route to the generation of surface waves and to outline an analytical model 
for the TSWD. Therefore, for the sake of simplicity, we use a cold fluid, non-relativistic 
plasma model that can be tackled analytically. However, we believe that the TSWD is a 
general process that may occur for laser and plasma parameters beyond the limits of the 
approximations used in the present paper. Relativistic and kinetic effects, the inclusion 
of damping mechanisms and the fully nonlinear evolution of the TSWD are left for future 
investigations. 

II. THE MODEL 

In our model, we consider a cold plasma with immobile ions and a step-like density 
profile Tii(x) = riiQ(x). The laser pulse has frequency u, linear polarization and wavevector 
k = (uj/c) (x cos 9 + y sin#). Translational invariance along z is assumed. For both s- and 
p-polarizations, all fields oscillating at the frequency u can be determined inside and outside 
the plasma from general Fresnel formulae for refraction and transmission at a boundary 
12|| between vacuum and a medium with an (imaginary) refractive index n given by n 2 = 



1 — u 2 /u 2 < 0, Up = (Annie 2 / me) 1 1 2 being the plasma frequency. 

In what follows we study the excitation of "H" ESWs with the magnetic field in the 
z direction. Thus, we deal with an effectively 2D geometry, in which B = Bz and the 
Maxwell-Euler equations read 
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V-E = 47re(n -n e ), (4) 

VxE = --d t B z , (5) 
c 

47T 1 

V x B = — j + -dtE, (6) 

c c 

m e d t v = — m e v ■ Vv — e + — x (7) 

and the current density is given by j = — en e v. 
We adopt the following expansion for all fields 

F(x,y,t) = F { (x) + eF {x,t-y sin 6/c) + e 2 [/+ (x, y, t) + /_ (x, y , £)] , (8) 

where e is a small expansion parameter and / stands either for the electron density or the 
velocity or for the EM fields in the (x, y) plane. In this expansion, Fi represents unperturbed 
fields of zero order (e.g., ni). The term F represents the "pump" field at the frequency u , 
that is taken to be of order e and is written as 

In what follows the cases with uj =uj and Uq = 2u are studied separately. The last term in 
(§) is the sum of two counterpropagating surface modes, which are assumed to be of order 
e 2 and are written as 

f±= l -[ft ± \x)e^y- i ^ t ]+c. C .. (10) 

Using expansion (||), the coupling between the pump and the surface modes is of order e 3 . 
Thus, from eqs.(f|-(5|-|6|-|7]) we obtain to order e the pump fields. To order e 2 , we obtain a 
set of linearized, homogeneous 2D Maxwell-Euler equations which have ESWs as normal 
modes. Note that terms such as v • Vv in the Euler equation are also of order e 2 , but are 
non resonant with the ESWs at the frequencies uj±. Such terms represent a source term for 
harmonic oscillations at 2uj and can be neglected in our model. 

ESWs are "H" waves which can propagate in a dielectric medium along a layer of dis- 
continuity of the dielectric function, with the latter changing sign across the boundary [13 . 



For a cold, step-boundary plasma there is no volume charge density perturbation associated 
with SWs, i.e. V • E = — AneSrig = 0. At a vacuum-plasma interface, electrons do not enter 
the vacuum side, but are "stopped" at the interface (x = 0) forming a surface charge layer. 
The dispersion relation for an SW of frequency u s and wavevector k s is 



u2 = uj 2 - uj = u 2 a s - 1 



(11) 



We have set for convenience a s = uj 2 /uj 2 . The dispersion relation ( |TT1) is shown in Fig.|T]. 
The EM field envelopes (see eq. fllOD ) are given by 



E s y (x) = E s y (0) [Q(-x)e q < x + e{x)e~ q > 3 
B s z {x) = l -^E s y (0) \®{-x)e q < x + e{x)e- q > 3 

E s x = ikE s y (0) Q(- x )—-e(x) 



-q > x 



q< 



q> 



(12) 



where g> = (u> s /c)(a s — lj\/a s — 2), q < = (u s /c)(l/ \Ja s — 2). The velocity fields are found 
from v s = (ie / 'm e uj S )E S . 

All feedback effects of the nonlinear coupling on the pump fields are neglected. The 
coupling terms that lead to the excitation of the ESW with frequency u± may be represented 
by the nonlinear force 



t+ L) (x)e ik±y - w±t 



m e (v T • Vv + v • Vv T ) + -(v x B T + v x B ; 



(13) 



In the term in square brackets on the r.h.s of ( |T3| ) only "resonant" terms with the same phase 

of 4 iVL) 

are included. Obviously, these terms exist if the matching conditions ([$]) hold. 
The nonlinear coupling leads to the parametric excitation and growth of ESWs. We thus 
let the ESW envelopes vary slowly in time, i.e. 



f±{x) -> f±(x,et) 



(14) 



To evaluate the growth rate, we use an energy approach. The surface energy density of the 
ESW that will be needed in the calculations is given by 
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Z7T J-iv/k J-oo x ' 

where the brackets denote average over one laser cycle and u± and are the volume 
densities of the kinetic and EM fields energies, respectively: 

^ = ^N', (16) 



M « = ±(|E ± | 2 + |B ± J 2 ). (17) 
The kinetic energy contribution vanishes for x < 0. Using eqs.(|l2"D and 5± x (0 + ) = 



v± >y (0) / y/a — 1, and finally summing the two contributions we obtain for the total surface 
energy U± = M±\v±, y \ 2 /2 where 

m e nica±(a± - 2) 1 / 2 (a| - 2a± + 2) 
M± = (^1? ' (18) 

To obtain an expression for the temporal variation of U± due to the nonlinear coupling 
forces, we write the Euler and the Ampere equation for the ESW fields keeping terms up to 
order e 3 and disregarding terms that are non-resonant with the oscillation: 

m e d t v ± = -eV ± + e^ L \ (19) 

V x B ± = — (-en oV± + eJ ( ± L) ) + -d t E ± , (20) 
c K ' c 

where is given by eq.(|l3|), and J± Vi ' ) = — e5n < f°^v±. Using eqs. (|l9| - p0|) together with 
Maxwell's equations we obtain for the temporal variation of the energy densities 

d t uf ] = -en D v ± • E± + eff L) • v±, (21) 
dtu^ = -V • S± - ejf L) ■ E±, (22) 

where S± = (c/4ir)'E± x B± is the Poynting vector of the ESW. We integrate the sum of 
the energy densities over space and average over time. Noting that the total flux of S± 
vanishes because of the evanescence of the ESW fields, we finally obtain the equation for 
the evolution of the total energy U of the ESW: 
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d t U± = — dy dxdt{——\v±\ + — ( E± + B± 

2,71 J-n/k JO \ 2 87T V 

= e- ^ rfx y ^ ^ rfy (v ± • (e5n^) E± + f^) ) . (23) 

The integral extends for x > only since there are no electrons for x < 0. 

From eq . (p3|) , performing the integral and neglecting high order terms one obtains two 
coupled equations in the general form 

^7^dt\v±\ 2 = A u\v + \\v + \S± sin if, (24) 

where A is the amplitude of the pump mode in dimensionless units, S± are positive factors 
depending on u p /u and 9, and if is obtained from the phase factors of the pump and ESW 
modes as <p = (p + + <p- — <fio- From (E3) one easily obtains 



8^\v ± \ = ±^ V \v ± \. (25) 

Setting sin (p — 1 simply corresponds to finding the relative phase of the growing modes with 
respect to the phase of the pump mode. The growth rate of ESWs is then 



S+S. 



T = A ujJ—± . (26) 

In the next section we use the analytical method outlined above to study TSWD in two 
cases of particular relevance: TSWD driven by the v x B force for normal laser incidence 
(which we name "2cj — > uj + uj" decay), and TSWD driven by the electric force for oblique 
laser incidence and p-polarization ("a; — > u/2 + u/2" decay). 



III. "2w -» uj + w" DECAY 

To discuss the u 2u — > u + u" decay, we restrict ourselves for simplicity to the case of 
normal incidence, as was done in the simulations reported in Ref. |[L0|| . Thus, the ESW 
frequency u± = u and k + = —k- The case of oblique incidence is indeed very similar, the 
main difference being that the matching condition k + + k_ = (2u/c) sin# holds, causing a 
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shift Sou 7^ between the two ESWs. It is found that the growth rate decreases monotonically 
with increasing incidence angle, since the pump force has a maximum at normal incidence 

0- 

The laser wave can be represented by a single component of the vector potential, A z = 
A z (x, t), that for x > is given by 

A z (x,t) = A z (x)cosut = A z (0)e- x/la coswt, (27) 

where l s = c/(u 2 - u 2 ) 1 ' 2 is the screening length. Imposing boundary conditions for the 
incident and reflected waves one finds 4,(0) = 2Ai(ujl s /c)(l + u 2 l 2 Jc 2 )~ 1 / 2 , where A { is the 
amplitude of the incident field. 

Electrons perform their quiver motion in the z direction. Thus, the magnetic force term 
is in the x direction and has a secular term (Ooj), named the ponderomotive force, and an 
oscillating term (2u) that, in what follows, we simply name the v x B force. The secular 
term corresponds to radiation pressure and creates a surface polarization of the plasma. 
The v x B force drives a longitudinal, electrostatic oscillation that acts as a pump mode 
for TSWD. In the expansion (^), the pump amplitude is supposed to be of order e. As will 
be shown below, coupling between ID and 2D fields occurs only in the overdense plasma 
(x > 0). Thus, since the v x B force is quadratic in the laser field, the expansion (H) 
also implies a(0) ~ e 1 / 2 , where a(0) is the (dimensionless) laser amplitude at the surface 
of the plasma, e.g a(0) = (eA z (0)/mc 2 ). For overdense plasmas a(0) ~ (u>/uj p )ai < a i: 
where a, = (eAi/mc 2 ). This yields e ~ [uj/u p ) 2 a 2 = (n c /n e )a 2 , where n c = m e uj 2 /Aiie 
is the "critical" density. We therefore expect our expansion procedure to be valid even at 
relativistic fields amplitudes a« > 1 for high enough plasma densities. 

We now derive the solutions for the ID pump fields at 2uj. The electron oscillation 
velocity is obtained from the conservation of canonical momentum as m e v z = eA z /c. For 
normal laser incidence, the longitudinal v x B force is given by 

2 

- -v z B y = e —A z d x A z = F°{x){l + cos2cjt), (28) 

c m e c z 

9 



where we have set F°(x) = {rrieC 2 /2l s )a 2 e 2x l ls = F°e 2x l la . To order e, one obtains the 
following equations for the longitudinal, electrostatic motion : 

rriedtV^ = -e{Ef ] + E {2w) ) + F°(x)(l + cos2cut), (29) 
dM 0) + 4 2UJ) ) = -4*e(6n<p + 8nf^), (30) 
d t 5nM = -nod x V^\ (31) 

All fields in the equations above decay inside the plasma as exp(— 2x/l s ). The secular part 
simply gives eE^(x) = F°(x) and 5n^ = —d x F°(x)/4:7re 2 = F°(x)/2ire 2 l s . For the motion 
at 2a; one obtains 

T>(2^) _ ~ tF ° ( V_ \ p -2x/l s ( on\ 

K -com e D[co 2 ) e ' {62) 
6h^= ni J° e~ 2x / l % (33) 

eE^ = ^-e~ 2x/h . (34) 
2D 

The denominator D is given by 

D=l -% = 1 -^, (35) 

which shows the well-known resonance at n e = 4n c due to excitation of plasmons with 
frequency u p = 2uj by the v x B force. 

The difference between the total number of electrons and ions for x > is, to order e, 

F° / cos2wr 



AA/" (l>0) 



/> K »» + , n f)^( 1 + ^). (36) 



The fact that AN^ X>0 ' > during most of the oscillation implies that, due to compression 
from the ponderomotive and v x B forces, electrons leave behind a charge depletion layer 
of thickness £ = AN^ X>0 ^ / 'rij. We note that £ is of order e, thus to lowest order it is correct 
to treat the charge depletion layer as a surface layer. Heuristically, the oscillating behavior 
of £ describes the plasma "moving mirror" effect ||15|| : due to charge depletion the laser 
is reflected at x = ( rather than exactly at x = 0. This leads to the appearance of high 
harmonics in the reflected light. 
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Since \D\ < 1, there is always a time interval in which electrons are pulled into vacuum 
forming a cloud of negative charge. This interval is very short for n e /n c ^> 4, i.e. D ~ 1, 
but for lower densities it is relevant in the interaction process. The motion in vacuum 
is strongly anharmonic and more difficult to solve analitically than the motion inside the 
plasma. However, it is found in section |TTT] that the surface modes gain energy during the 
phase of electron motion inside the plasma only, so that the expressions of fields for x < 
are not needed. Nevertheless, we have to assume that the density of the electron cloud for 
x < is low enough for the laser propagation not to be affected. 

Using eq.(32), the NL force ( |i~3"D may be written in terms of the oscillation velocity V^ 2 ^ 
as follows: 



cNL 



m e (y^\x)d x v T + v x ^d x V^\x)*) - -V^(x)B ZjT y 



(37) 



Because of the symmetry with respect to inversion of the y axis, the two ESWs may 
differ only by a phase factor and have the same amplitude and energy. Inserting the force 
( P7|) in the energy equation (^) and performing the integral in y we obtain for the ESW 
energy U = U± 



d t U = \ E r°° dx Ki ■ (e$n e ^ (x)E* , - m e n v* x .^V^ (x) 



l=+k,-k' 

-m e n V x ^ ) (x)d x r i + n -V x ^\x)B* z _,y ) + c.c. 

c ' 



(38) 



We rewrite all fields as functions of V^ 2 ^ and v±. We obtain, after some algebra 



m e rn 



-m. I vS 2u3) 



2{l + q> l 

All the terms of eq.(|39D are proportional to 



1 + q> l s )v* + v*_ + 2v*v* 



(39) 



3? (y( 2uj) v* + v*_) = \V x {2uj) v 2 ± \ cos(0 + tt/2), 



(40) 



where = (f> + + 0—, <j)± being the phase angles of v±, and we used eq.(^2]) and the fact that 
F° is real and positive. Thus, the phase of the growing ESWs with respect to the v x B force 
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is such that <fi = — 7r/2, the value for which the growth rate is positive and has a maximum. 
The overlap of the two ESWs gives 

^+ iky—iu)t I ^~ —iky—iujt l^i I —iujt f iky+i<f>+ _i_ —iky+i<f>-\ 

2 2 ~ 2 ^ / 

= |C ± |e-^- i7r/4 cos (ky + A0) , (41) 

where A<p = (<ft + — 0_)/2. This is a standing wave which has a temporal phase shift — 7r/4 
with respect to the v x B force of eq. (|28|) . The phase is such that, at a given position in y, 
the temporal maxima of V^ 2 ^ and v ^ overlap once for laser cycle. The angle A0 gives the 
location of the maxima of the standing wave on the y axis, which depends on the arbitrary 
choice of the initial phase. 

Using eq. ( ^0|) , we can rewrite the energy variation as 



d t U = m eni \V^ 



in) I 



l 2 + 



i + q>l 3 q>l s (i + q>h 



(42) 



Note that d t U > since g>/ s = y (a — l)(a — 2) _1 > 1, where a = oj^/uj 2 = n e /n c . 
Eliminating V^ 2 ^ as a function of a, and a, after some algebra the growth rate T is obtained: 



2 



[a 



- If' 2 



o;|a-4|[(a-l) 2 + l](a-2) 1 / 



2 



2(a- 1) 



a(l + q>l s ) ag>/ s (l + g>/ s 



- 4wa, 2 — 777-^ -ryr. (43) 

4 a|a-4|[(a-l) 2 + l](a-2)V2 v 7 

The leading contribution was highlighted in the last equality. 

A plot of the growth rate as a function of a is given in Fig.[| For moderately overdense 
plasmas, T is a considerable fraction of the laser frequency. The resonance at n e = 4n c 
is due to the resonant behavior of the pump field. The growth rate diverges also in the 
electrostatic limit when n e —>■ 2n c . However, in this limit the value of k tends to infinity, 
i.e. the wavelength becomes very small and one expects that this second resonance might 
be damped by thermal effects, which are neglected in our model and are left for future 
inevstigations. 

The comparison of spatial and temporal scales predicted from this model with the sim- 
ulation results shows reasonable agreement for laser and plasma parameters such that our 
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expansion procedure is valid |T0| . We note that simulation results for high intensities in the 
relativistic regime suggest that a 2uj —>■ uj + uj TSWD-like process still occurs and produces 
strong rippling of the plasma surface; the spatial scales are different from those predicted 
by our analytical, cold fluid model, but still of the same order of magnitude. 

We conclude this section by noting that, although we took the laser polarization along 
z, the whole derivation is valid at normal incidence for any polarization direction. The only 
difference in 2D geometry is that, for the laser polarization along y, the quiver oscillations 
at uj and the ESW oscillations at the same frequency overlap. This is observed in numerical 



simulations [14 



IV. "u ^u/2 + u/2" DECAY 

A case of particular interest is that of oblique incidence and p-polarization of the laser 
pulse, since these conditions lead in general to a stronger plasma coupling. In the preceding 
section we have found that the 2uj — > uj+uj decay is driven by the longitudinal component of 
the velocity field, as shown by eq. (p7|) . For p-polarization and 8^0 the electric field drives a 
strong longitudinal oscillation at the frequency uj; thus, one expects that the TSWD process 
leads to the generation of subharmonic ESWs, having frequencies around uj/2, with peak 
efficiency at some angle of incidence. As shown below, this picture is substantially correct. 
However, there is a small but non- vanishing growth rate for the uj — > uj/2+uj/2 process even 
at normal incidence. In fact, the pump force (|i~3| ) is now given by 

ff L) = -m e (v^d x VM + V^d x v T ) - me^VM + V y ^d y v T ) 

-^(V (w) x Bzp + Vzp x B (w) ), (44) 

and does not vanish for 8 = 0. 

The case of oblique incidence and p-polarization of the laser pulse may be tackled with 
an approach analogous to that of section |TT1[ The main physical difference with the previous 
case is that now the "pump" is a divergence-free velocity field at the laser frequency u, and 
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its amplitude is now proportional to the laser field rather than to the laser intensity. In 
applying the expansion procedure (§) again, now e ~ {uj/uj p )ai (appendix [A]). Thus, the 
limits of validity of this expansion are more restrictive in the present case: the expansion 
tends to be valid only for rather high densities or rather low fields. One must also note 
that temperature effects, neglected in the cold plasma approximation, are important for low 
fields. 

The calculation of the growth rate for subharmonic ESWs can be performed by the same 
method as in the preceding sections in a straightforward way. The details are reported in 
appendix |A[ As a final result, the growth rate V = T(a, 9) may be written as 



where Tb is the Fresnel factor for the magnetic field and the factor K(a, 9) scales weakly 
with density. A plot of T as a function of 9 for different values of a is reported in Fig||. 
The frequency shift 5uj of the two ESWs can be calculated as a function of 9 and a from 



Fig.|3], one finds that 5u ~ O.lu in conditions favorable to the uo — > uo/2 + uj/2 TSWD. 



The TSWD concept has been investigated analytically in two cases that clarify its basic 
features in the context of the cold fluid plasma approximation. This process may lead to the 
excitation of ESWs by an ultrashort laser pulse in a solid target without any special (e.g., 
grating-like) structure, and even for normal incidence and s-polarization. 

In both the cases discussed in this paper, we found a strong decrease of the TSWD growth 
rate for increasing densities. This does not necessarily imply that TSWD is not relevant to 
short pulse interaction with solid targets which have n e /n c ^> 1. In fact, important processes 
such as high harmonic generation or fast electron production are more efficient when the laser 
pulse interacts with a moderate density plasma. This is the case for most of the experiments 




(45) 



the matching condition k + + k_ = (uj/c) sin0. The result is shown in Fig.[|. Comparing with 



V. DISCUSSION 
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of short pulse interaction with solid targets, since a moderate density "shelf" is usually 
produced at the time of peak pulse intensity by target ablation and plasma hydrodynamic 
expansion during the leading edge of the short pulse or during the long prepulse. For 
instance, the ll 2u — ► uj + uj" decay appears to be very efficient exactly in conditions favorable 
for high harmonics production via the "moving mirror" effect [|T3|,|TBJ. Simulations fTU| 
suggest that TSWD in strongly nonlinear and relativistic regimes (beyond the limits of the 
analytical approach of this paper) may produce surface perturbations acting as a "seed" for, 
e. g. , electron instabilities leading to current filamentation |l7|Jl8| or detrimental distortions 
of the moving mirrors, and for Rayleigh- Taylor-like (RT) hydrodynamic instabilities occuring 
on the time scale of ion motion. In Ref. |H| it is estimated that the growth rate of the RT 
instability is much slower than that of TSWD. A different hydrodynamic instability with a 
surprisingly high growth rate has been reported in Ref. fL9f). 

The 2uo — > uo + uo and the uo — > uo /2 + uo /2 processes have been considered independently. 
This is appropriate since these are resonant processes which do not interfere with each other. 
In principle, both decays may occur during the interaction of an intense laser pulse with 
an overdense plasma. According to our analysis the 2uo — > uo + uo process appears to have 
a stronger growth rate. However, in the case of the uo —> uo/2 + u/2 process our analytical 
approach is valid for a very narrow range of parameters only. Either an extension of the 
present analytical approach or numerical simulations are needed to investigate the TSWD 
for a very intense, p-polarized laser pulse at oblique incidence. 



VI. CONCLUSIONS 

We have discussed the Two-Surface Wave Decay as a novel mechanism for the excitation 
of electron surface waves in the interaction of ultrashort, intense laser pulses with overdense 
plasmas. TSWD is based on the parametric excitation of a pair of ESWs. The "pump" force 
may either be the magnetic or the electric force of the laser pulse, and leads respectively 
to the generation of two ESWs with frequencies close to the laser frequency of half of it. 
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An analytical model for TSWD in the cold fluid plasma, non-relativistic approximation has 
been developed. The model supports the interpretation of recent simulation results [TIJ and 



suggests that TSWD may be of relevance in certain regimes of laser interaction with solid 
targets. 
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APPENDIX A: GROWTH RATE OF THE uj -> uj/2 + u/2 DECAY 

We now give the detailed derivation of the growth rate for the u — > u/2 + u/2 TSWD of 
section [TV} The pump fields are easily found with the help of Fresnel formulae and Maxwell 
equations. For instance, the magnetic field is given by 

B^{x,t) = B ^\0 + )e ikys ' me - x/lp - lujt + c.c. , (Al) 

where the screening length is given by 

*P = - ( 1 - -2 cos2 e ) V2 = ~ , 1 2a > (A2) 
w p \ / uj Va - cos 2 9 

and the magnetic field at the surface is given by the Fresnel formula 

B z (0+) _ 2n 2 cos# 
B z>i Vn 2 — sin 2 9 + n 2 cos 9 

= 2(a ~ 1)C " Sg =F B {0), (A3) 

(a — 1) cos 9 — i\f a — cos 2 9 

where B z i is the incident field amplitude in vacuum. The electric field components are found 
by using V • E (uj) = and V x E H = i(u/c)B^ + c.c. 

The nonlinear force is given by eq. (p7|) . Evaluating the spatial derivatives as d x \^ = 
-V^/I p , dyV^ = iktV, d y v± = ik±v±, and d x v± = —q±v± where q± = q > (u±), and 
keeping resonant terms only we find 
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r(NL) _ 
J±,x — . 



1 



f: 



(NL) _ rne 
±,v ~ a 



12 J ff~\j gC 

1 ~, 
v 

i(k t - k T )V^v^ y 



( V {uj) Bt + B { 



z u t,v 



e -(q T + l/l v )x + QQ 



(A4) 



{q T +l/l v )x 



] B* + B^v* 

x T,z ' z =F,a; 



+ c.c. (A5) 

We rewrite all "pump" terms as a function of V^ 03 ' by using t/. = ik t l p V y ^\ B^ 3 ' = 
(m e c/el p )(k^lp — 1)V1 ■ in a second step we rewrite all SW fields as a function of v± iV by 
using v± )X = (ik±/q±)v±^ y and b± )Z = —m e uj±(a± — l)v± ty /(ecq±). We thus obtain 

,2 



k lo 1 I 

-± (A* + fc^Zp + fc T ) + ^(« T - 1) + y 



— a v y u T,y^±> x(i 



7{NL) _ me™. t 



=F,1T 



, , / 1 hk T 
hip \ q T ~ y 



-{q T + l/l p )x 



(a T -l))+k T 



a v y D ^,y l ^±,y (i 



The relations u>±(a±—l) / q± = u±c^/a± — 2, k±/q± = 1/ \Ja± — 1, k t l p = sin6 l (a— cos 2 
may be used to simplify the expressions for Q± >x and Q± >y . We thus obtain 



Q 



±,X 



V 7 ^ 



-\k t {\ + q T l p ) + k T ] + ~-^Ja T - 2 + 



1 



Q±, y = hip [q T - j- 



^ x la T -2) +k T . 



(A6) 
(A7) 
(A8) 
(A9) 

0)-l/2 

(A10) 
(All) 



\ l P V"t 

In the energy equation (p3|) , since (faA^ = the integrand, eliminating w± iX , is given by 



v ± ■ £') = ^V^vlA, (^Q ± „ + Q ± ,„) + ex. 



Q± 



~-Q±,x + Q±, 



We thus find 



d t U± = n. 



-Q±,x + Q±, 



simp. 



(A12) 
(A13) 

(A14) 



1 + (q+ + q-)l p \q± 

Here, <p = (<j> + + </>_ — $), where $ and 0± are the phase factors of V y ^ and v±, y , respectively 
It is worth rewriting these equations in a more compact form. Introducing eij = eBi/m e ujc 
and using Fresnel's formulas to eliminate V y we obtain 
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Tfl 

d t U± = ni^a i c\J r B {9)\\v+,y\\v-,y\suiipG±(a,0), (A15) 



8 

where we have posed G±(a, 9) = G (a, 9)g±(a, 9) and 



Va - cos 2 9 

°° ~ (a-l)[l + ( g+ + g _)Z p ]' (A16) 



9± — lp \ —Q±,x + Q±.y 



a+-l)(a_-l) c V"± " 1 V"± " 1 



+ktl p \q^l p - 1 - ^=%y " V a T " 2J + Mp- (A17) 

Using the equations above we obtain the following coupled equations in the form (|4]) for 
the field amplitudes 

,y\ — ai^\v+y\\v^y\\J-B\G± sin ip, (A18) 

where fi± = 4M±/(m e nic). Thus, the growing modes have relative phases such that simp = 
1. The growth rate is given by 



T = aMF B \G > = aiw\r B \ ( —) K(a, 9) (A19) 

\ \n e J 

which is the formula (f£5|). A plot of T is depicted in Fig|| 
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FIGURES 

FIG. 1. Dispersion relation of electron surface waves (thick lines) eq.(^TJ), and the matching 
conditions for TSWD, eq.(||). 

FIG. 2. The TSWD growth rate T (normalized to af ui) for the 2ui — > ui + ui process at normal 
incidence, eq. (|43|) , as a function of a = n e /n c . 

FIG. 3. The TSWD growth rate T (normalized to aiui) for the u —* uj/2 + uj/2 process at 
oblique incidence, eq. (fj5[) , as a function of 6 and a = n e /n c . 

FIG. 4. The frequency shift 5uj for the uj — > to/2 + uj/2 process at oblique incidence, as a 
function of 9 and a = n e /n c . The dashed line gives the shift 8oj max = (1/2) sin 9 that is obtained 
for a — > oo. 
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